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ABSTRACT 

A generic physical situation is considered where Im 11, the imaginary part of polarization 
operator (generalized susceptibility), can be measured on a finite interval and the high 
frequency asymptotics (up to a few orders) of 11 can be calculated theoretically. In such 
a case, it is desirable to derive an equivalent form of the Kramers-Kronig dispersion 
relation, the so-called sum rule, in which both the high-frequency part of Im 11 in the 
dispersion integral and the high-order contribution to 11 are suppressed. We provide a 
general framework for derivation of such sum rules, without any recourse to an infinite- 
order differential operator. We derive sum rules for a wide set of weight functions and 
show that any departure from the e~* behaviour of the weight function in sum rules leads 
to modulating factors on the theoretical side of sum rules, providing its low frequency 
regularization. We argue that by including modulating factors one can extend the domain 
of validity of sum rules further to an intermediate region of frequencies and can account 
for "bumps" which were observed numerically on the phenomenological side of sum rules 
at "intermediate" frequencies. 
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1. Introduction - Kramers-Kronig dispersion relations for generalized susceptibility 11, 
such as in 

U{tu}) = - / — — - sds, (1) 

IT Jo S'' + 

where a; is a frequency, arise in many branches of physics. The reason is that disper- 
sion relations are direct consequence of the analyticity of generalized susceptibility in the 
upper-half complex plane which, in turn, is believed to be the consequence of causality. 

Generalized susceptibility 11 determines the response of the system to a perturbation. 
Dispersion relation (|1]) relates the value of 11 at imaginary frequency to the integral over 
its imaginary part, Im 11. In the generic physical situation to be considered here, Im 11 is 
determined experimentally on a finite interval (0, c^o) via directly measurable quantities 
such as scattering cross sections. On the other hand, we shall suppose that a few orders 
of the asymptotic expansion of 11 for u ^ 1 can be calculated theoretically. An example 
is provided by quantum chromo dynamics (QCD), where the generalized susceptibility is 
the polarization operator, defined via the current-current correlation function, 

n^^(g) ^ i J e-'^^{0\ Tinx)fm |0) d^x = {g^^^ - q^q^ n(g2), (2) 

and g is a four-momentum. The polarization operator satisfies the standard QCD disper- 
sion relation 0] 

MQ') = - , , XL ds, (3) 



where Il{Q^) = n(— g^). The imaginary part, Imll, is proportional to measurable cross 
sections such as that for e~e+ (electron-positron) annihilation into hadrons and, for suffi- 
ciently small momenta, can be supplied from experiment. On contrary, one can get, for 
sufficiently large, the QCD representation for the polarization operator as (asymptotic) 
series. 

Here, the first term oq corresponds to asymptotic freedom and the remaining coefficients 
an {n > 2) parametrize power corrections, describing long-distance effects. They are 
given as the product of a short-distance term calculable within perturbation theory and a 
corresponding (quark, gluon, etc.) condensate which is unknown in theory 0. Therefore, 
in the generic physical situation, it is desirable to find an equivalent form of (|T]), a sum 



rule, in which both the high-energy part of Im 11 and the high-order contribution to 11 are 
suppressed, thereby reducing respective contributions of Im 11 and 11 for those frequencies 
and those orders for which our knowledge is incomplete. 

The first step in this direction was performed by Shifman, Vainshtein, and Zakharov 
(SVZ hereafter) [0]. By applying the formal differential operator Lm, 

(;r^(«0"(-5^)""'. w 
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on both sides of dispersion relation @ they derived the following equivalent form of the 
dispersion relation: 

if e-/«Vn(.).,^g^i:-^, (6) 

called the SVZ sum rule. In the contribution of Imll to the dispersion integral is 
reduced for energies s > Q"^ hj an exponentially decreasing factor. On the right-hand side 
of the contribution of high orders of perturbation theory for 11 is suppressed due to the 
factor n\ in the denominator. Therefore, using sum rule (|^), theoretical predictions and 
experimental results can be tested in a more efficient way than using original dispersion 
relation (|I]). In particular, in the case of QCD, one hopes that the validity of sum rule 
d^) can be extended to intermediate energies and that the condensate values could be 
determined. A natural question arises whether one can provide a general framework for 
derivation of sum rules with a general weight function and specify their domain of validity. 

2. Results - To answer the above questions, we considered the class of weight functions 
of the form 

= (l/«)t(^/")-^e-*^^^ (7) 

where < a < P, and 

X,{t) = t^e-^' (8) 

with 7 > 0. Our choice of weight functions covers a wide range starting from 

monotonically decreasing to "Gaussian-like" (cf. Ref. [§) weight functions. 

Our main result (see the next section) is that a general sum rule with a weight function 
x{t) can be written as 



- rx^'\s/Q')lraTl{s)ds^Q'Y.{-^r''\^''\eJQ')^^ 
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where x''"'^ denotes the nth derivative of xi^)- Constants £„ in @, < ^ Q^, are, in 
general, different for different n. 

We shall show that physically reasonable sum rules arise only if x{t) together with all 
its derivatives is regular at the origin. Then, in the region <^ Q^, 

^ X^'\smimUis)dsr~.Q'J2i-l)(--^^Xn^, (10) 
vr Jo 

where the coefficients Xn are defined by 

oo 

x(t) = ExX. (11) 

n=0 

Sum rule ([T0| ) includes all particular cases discussed previously in Refs. 0, ^. Because 
weight functions are entire functions, the coefficients Xn as a function of n decreases to 
zero faster than any polynomial in the limit n — > oo. 

Sum rule (|^) will be derived by employing theory of summability methods, without 
any recourse to infinite-order differential operator such as (^. This approach will also 
enable us to include into our consideration the important case of weight function Xj- In 
contrast to e~* where frequency is required to be imaginary, the weight function x^y allows 
one to consider dispersion relation (|l|) for any frequency in the upper-half complex plane. 

In Sec. 4, we shall show that unless x{t) = one has to have e„ > 0. Thus any 
departure from the e~* behaviour of the weight function in the sum rules leads necessarily 
to Q^-dependent factors X^'^K^n/Q'^) on the theoretical side of the sum rules, which we 
shall call modulating factors. Sum rules as a function of parameters a, (3, and 7 [see 
Eqs. (^ and (H)] for a particular weight function are analyzed in Sec. 5. Domain of 
applicability of sum rules (^ and (P) and role of modulating factors are then discussed in 
Sec. 6. 

3. Derivation of general sum rules - In what follows, the term "summability method" 
will stand for a moment constant analytic (Borel-like) summability method 0, ^- In 
the theory of summability methods, weight functions Xo,/3 and x-y are examples of the 
generating function x{t) which generates the moments n{n), 

li{n)= / xit)t-dt. (12) 
Jo 

Any summability method is characterized by its generating function which, obviously, 
must decrease to zero faster than any power of t in the limit t ^ 00. The properties 
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of a summability method then depend on the properties of the corresponding moment 
function F{t), determined in terms of moments fi{n) 0, 

oo ±12 

In the case of Xa,/3, the moments are /i(n) = T{an + 13) and the corresponding moment 
function is known to be the Mittag-LefHer function 

oo n 

In the case of X7, the moments ^{n) are growing roughly like (In ra)" and the corresponding 
moment function is analyzed in [Q. Throughout the paper, we shall speak of Borel 
summability if and only if x(t) = e~*, in which case the moments are ^{n) = n\. All other 
choices of x(t) will be referred to as the other summability methods. 

For any weight function x{t) discussed here, the Cauchy kernel in dispersion integral 
(1^) can be expressed as 

roc roc 

/ x{Q't)F{-st)dt= x{st)F{-QH)dt, (15) 
Jo Jo 



where integrals in (^) are absolutely convergent P, ^. Representation ([TsD is the key 
relation which will allow us to represent the Cauchy kernel in dispersion integral (|]) as a 
X-weighted integral and, at the same time, to find the coresponding modification of the 
right-hand side of (^. 

In the next, we shall assume dispersion integral (|]) to be absolutely convergent [§. 
Otherwise, one can always consider, instead of (^, higher derivatives of the dispersion 
relation, 

m') = - I T—^^ds. (16) 



n\ \ dQ^ J ' 7lJo{s + g2)n+l 

Relations of the type (|T6|) do not bring any complications to our discussion, since absolute 
convergence of the integral in ([T^) allows us to differentiate within the sign of integration 
and represent any power of the Cauchy kernel as 

1 1 I d\ 1 1 f'^ „, ^r, , d 



(s + g2)n (n-1)! \ ds s + g2 {n 



(17) 
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Again, for any n the integral in ([T7| ) converges absolutely Then, by using relation 

(|T6| ) for = 1, dispersion relation can be written in the form 

Now, the key problem is to find the function G{l/M'^) such that the right-hand side of 
([T8|) can be represented as an integral of the form 

E 0^ - r F(-QVM-)G{llKf-) . (19) 



n=0 



In order to find (^(l/M^), one again uses representation (|T5D and (|T7D of the Cauchy 
kernel and its powers. Since is assumed to be sufficiently large, one introduces only 
small error if is approximated by 1/ [Q"^ + e)", where < £ <^ Q^. Then, according 

to (ll^), general power can be approximated as 

^n— 1 

Q2n [Q2+sY (n-"T)!7o 

which introduces an infinitesimally small correction to higher orders. In what follows, we 
shall first approximate the right-hand side of (|18]) according to (^O]) and we shall discuss 
the limit £ afterwards. Now, in the ^-approximation, 

G(l/M2) = E(-l)"x("n^/M^)n^, (21) 



roc 

-^j^ F{-QH)x^-''\et)r-'dt, (20) 



and comparison with Eqs. (|T8]) and ( |T9| ) leads immediately to general sum rule (^. 

4. Modulating factors.- A natural question arises what is the role of these yet unspec- 
ified parameters e„. Can one get rid of them? 

In the case of the Borel method, the generating function x{t) and the moment function 
F{—t) are identical, 

xit)^Fi-t). (22) 

We shall refer to this property as self- duality" . It means that both x(t) and F{—t) are 
exponentially decreasing and integrals in (|1^), (0), and ( pOD are absolutely convergent 
Laplace integrals. In this case, obviously, e can be sent to zero in (pO|), 

G{l/M') = V (23) 
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and one recovers the SVZ sum rule (^. 

For all other moment constant summability methods discussed here, the "self-duality" 



(11) is lost, 

Xit) + Fi-t). (24) 



Even more significant than the lack of "self-duality" property ([2^ ) for the other summabil- 
ity methods is that, in the case of the Borel method, F{—t) (= e~*) approaches zero faster 
than the inverse of any polynomial in the limit t —>■ oo, while such inverse polynomial 
decrease is charactersitic for all other summability methods. Indeed, for all other summa- 
bility methods, the asymptotic behaviour of F{—t) in the limit t — > oo is characterized by 
an abrupt switch off from the exponential decrease e~* to a universal polynomial decrease 

i, 

F{-t) ~ 0{t-^) {t ^ oo). (25) 

For example, in the case of Xa,f3{t) the moment function is Ea^^it) [see (^^], having the 
asymptotic behaviour p 

oo / j.\n 

^1 r(/? - an) 

Unless both j3 and a are integers and (3 < a, polynomial terms in the asymptotic ex- 
pansion ( p6D of Ea^—t) are always present. Indeed, Eq. ( p6D implies that Ea^—t) is 
exponentially decreasing if and only if for all ra, /3 + an is a pole of the Euler gamma 
function T{z). For x-yif) see [^. 

The loss of the exponential decrease of F{—t) for all other summability methods implies 
that the limit e ^ cannot be taken. Because of the universal behaviour ( ]25| ) of F{—t) 
at infinity, 

I n-QVM^)ijp) d-, (27) 

diverges for any n > and G(l/M^) cannot be pure polynomial. It must contain modu- 
lating factors which ensure convergence of (|T^) for small M^. Therefore, in general case, 
the sum rule corresponding to the weight (generating) function x{t) must be given by 
relation (|) with < e„ < 

5. Analysis of sum rules for different weight functions - Let us first consider Xa,i3 and 
denote p = jS/a — 1 and q = 1/a. Then, in order to have meaningful sum rules, p and 
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q must be respectively nonnegative and positive integers. Otherwise, ~ Xa,(}{t) 

in the hmit t — 0, resulting in the fact that all orders in (|^) will behave as 0{1) in the 
limit ^ cxD on the theoretical side of sum rule and thus, they will merge to a single 
order and become equivalent. 

If p is a nonnegative integer and g is a positive integer, then Xa,pi including all 
its derivatives, is regular at the origin. Therefore, for e„ ^ Q^, modulating factors 
xi"^(^n/<5^) in (|^) can be approximated by their values at the origin and one finds rela- 
tion (|l^) which, upon appropriate choice of the weight function, gives the SVZ sum rule 
d^) and its generalizations discussed in |Q. 

To discuss particular cases, one uses 

;^W(0) = n\xn = (-l)'g ^^^5n,p+., = ("l)'? jKv+H- (28) 
Therefore, unless n = p + Iq, where / is a nonnegative integer, xi"fl(0) = 0. Relation ( pH]) 



implies that the sum rule with such weight function is (i) selective: it only selects the 
orders n = p + Iq of the asymptotic expansion (|p, (ii) the contribution of the {p + /g)th 
order in is then reduced by the factor l/l\ on the theoretical side of the sum rule. 
For example, in the case of x(^) = 2e~*^ only even terms of the asymptotic expansion 
(0) enter the sum rule (p!o|). Thus, such selective sum rules can be particularly useful in 
the case when some symmetries of the problem ensure that only {p + /g)th orders of the 
asymptotic expansion are nonzero. 

In the case of x-yy similar argument shows that it is useful to have 7 7^ only if the 
first 7 orders of the asymptotic expansion (^) are identically zero. Otherwise, for e~^*, 
one has 

Jn 00 / T\k 

-exp(-e*)^n!x„ = E4fr. (29) 

Estimates, using Euler-Maclaurin sum formula [0, give \xn\ ~ (lnn)"/n! in the limit 
n 00. 

6. Domain of validity of sum rules and physical significance of modulating factors- 
Sum rules (^) and (y) are derived for sufficiently large. However, the meaning of 
"sufficiently large" must be decided separately in each particular case. Since, eventually, 
one wants to make some practical use of sum rules, one is interested in some intermediate 
region of energies (or frequencies). The best possible situation occurs if there is some 
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nonempty intersection ("fiducial region") of respective intervals on which Imll and 11 
are known. Note that if the intersection is empty, one may gain by repeating all our 
calculations with x{t) replaced by XA{t) = x{^t)y by introducing a scaling variable A. 
This scaling leads to F{t) replaced by -Fa(^) = F'lAt) and n\ factor in (PT|) replaced by 
n\/A^. The corresponding sum rule is then 

1 /-oo 4" 

- / ;^a)(.A/Q2)n(.)ci.~Q2^(-l)(«-l)x(")(^nW)«n^7^- (30) 

71 Jo nlQ^"- 
For A > 1, one gains on the experimental side of sum rule: the interval on which Imll 
must be accurately determined is scaled down from interval (0, Q"^) onto interval (0, Q'^/A). 
However, one lost on the theoretical side of the sum rule: (i) damping of higher order terms 
is slower, (ii) since the factor A^/nl has its maximum at n = A, terms in the expansion 
(^) at the order n ^ A enter the sum rule with relatively the highest weight. Therefore, 
although the scale transformation is worth considering, it is not at all obvious whether 
one can extend the validity of sum rules to the intermediate region. 

We argue that modulating factors play important role in the extension of validity of 
sum rules to the intermediate region. The argument is as follows: by comparing the 
behaviour of both sides in (|^) in the limit ^ one finds that the left-hand side goes 
to zero while the right-hand side diverges. On the other hand, if modulating factors are 
present then, in contrast to (|^), both sides of (|^) tend to zero in the limit —>■ 0. 
Modulating factors can be viewed as the low energy regularization of sum rules. The 
mismatch between the vanishing rates of the two sides of (P) in this limit (the left-hand 
side tends to zero polynomially while the right-hand side tends to zero faster than any 
power of Q"^) is not too important: nobody expects sum rule (^) to be valid for too 
small. Nevertheless, the very fact that, in the presence of modulating factors, both sides 
of d^) tend to zero in the limit —>■ 0, represents substantial improvement over the SVZ 
sum rule. 

Another argument to support our claim that modulating factors can extend the validity 
of sum rules to the intermediate region is the appearance of "bumps" observed numerically 
on the phenomenological side of sum rules in |Q at ~ 0.5GeV^ (see Fig. |I]). In the 
absence of modulating factors, it was impossible to match the fits of both sides of sum 
rule (^) in the intermediate energy region and a conclusion was drawn that the Borel 
summability is optimal and that by any departure from the Borel summability sum rules 
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Q^EGeV^] 

Figure 1: Formation of a "bump" on the phenomenological side of sum rule for q = 1, 2, 
and 3. As q increases, the bump becomes more and more pronounced. 

get worse [Q. However, if modulating factors are included, then, since G{l/Q'^) ^ in 
the limit 0, G{l/Q'^) will develop a maximum at the point where X^'^~^K^n/Q'^) 

overcomes polynomial increase of l/Q^" (in the limit 0). This turning point then 

provides a way of determining modulating constant En phenomenologically. 

7. Discussion and conclusions - We have provided a general derivation of sum rules 
within the framework of summability methods, without any recourse to a formal differ- 
ential operator (D (cf. 0, 

We have shown that unless symmetries of a physical model cause some particular 
orders of the asymptotic expansion of the polarization operator 11 to vanish, the 
only reasonable choices of weight functions, from the set (|^) and we have considered, 
are e~* and e~^*. The use of e"*^* instead of e~* in sum rules means also a significant 
qualitative change: due to very special properties of the summability method based on 
the weight function e~^* [|ll[| , one is not bound to consider dispersion relation (|l]) for only 
purely imaginary frequencies as in the case of e~* but, instead, one can consider dispersion 
relation and derive sum rules for any frequency in the upper-half complex plane. This 
can have useful applications in other physical models. Moreover, a sum rule with e~^* has 
two advantages over sum rules with e~*: 

• the high-energy part of the integral is cut-off more effectively; 

• the onset of the cut-off in the integral (^) starts sooner, already at s ~ In 2 and 
not at s ~ as in the latter case. 
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All the above advantages may compensate for a slower, (lnri)"/n! damping [cf. Eq. p9| )] 
of higher orders of expansion (^) in comparison to l/n\ damping in the case of e~*. 

Further, we have proved that if a general weight function is used, the sum rule must 
contain modulating factors. The latter provide the low frequency regularization of sum 
rules and could extend their validity up to the intermediate region of energies. The 
very existence of the modulating factors in a sum rule with weight function different 
from e~* results from the universal 1/t decay [see Eq. of the moment function 

F{—t). It is amazing to note that a modulating factor can be also present in the case of 
the Borel summability method, and may lead to essential improvement of the fit. One 
cannot rule out its presence unless one has absolute control over numerical values of 
expansion constants in the asymptotic behaviour (1) of 11. If not, one can always find a 
positive constant e such that replacement of by !/(£: + Q"^) makes a change in the 
expansion (^) which is within the error due to our incomplete knowledge of the expansion 
constants a„. Then, the use of (pO|) inevitably leads to a modulating factor. The only 
exception of the Borel summability method is that, in principle, the modulating constants 
En can be sent to zero and the limit exists [giving the SVZ sum rule (U)], while for all 
other summability methods, the modulating constants must be nonzero. The modulating 
factors were not seen using the formal infinite-order differential operator Lm involving two 
limiting procedures [cf. Eq. ||)], because the latter set the argument of the modulating 
factors to the origin rendering them independent. 
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